Magnetic domains in III-V magnetic semiconductors 
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Recent progress in theoretical understanding of magnetic anisotropy and stiffness in III-V mag- 
netic semiconductors is exploited for predictions of magnetic domain characteristics and meth- 
ods of their tuning. We evaluate the width and the energy of domain walls as well as the 
period of stripe domains in perpendicular films. The computed stripe width d = 1.1 /im for 
Gao.957Mno.o43As/Ino.i6Gao.84As compares favorably to the experimental value 1.5 /jm, as deter- 
mined by Shono et al. [Appl. Phys. Lett. 77, 1363 (2000)]. 
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The utility foreseen for ferromagnetic semiconductors 
rests on the possibility of tailoring their electronic and 
magnetic properties on the same footing. Indeed, suc- 
cessful control of the Curie temperature T^-by hole 
density haj|-,been achieved in Mn-based IV-VlJil III-V,El 
and II-VHa semiconductor compounds. It has also 
been found that both biaxial strain imposed by lattice 
mismatchE and confinement in quantum structuresa can 
serve to engineer the direction of the easy axis. Further^, 
more, light-induced changes of the magnetic phase intll 
(In,Mn)As/(Al,Ga)Sb andS (Cd,Mn)Te/(Cd,Zn,Mg)Te 
heterostructures has been demonstrated. Finally, the 
possibility to tune Tq of (In,Mn)Aa quantum wells by 
metallic gates has also been shown.13 We can anticipate 
many new demonstrations of tunable magnetic and mag- 
netotransport properties in the future. 

These tuning capabilities along with Tq's as high asS 
110 K in (Ga,Mn)As with 5.3% Mn have triggered a con- 
siderable theoretical effort to elucidate the origin of fer- 
romagnetism in III-V magnetic semiconductors. While 
there is general agreement that the Mn constituent intro- 
duces both localized spins and itinerant holes, the nature 
and energy of the Mn-derived states, the role of intrinsic 
defects, the relative importance of charge and spin fluc- 
tuations, as well as the consequences of electrostatic and 
magnetic disorder are still under debate.tl 

Recently, a quantitative theory of hole-mediated ferro- 
magnetism in tetrahedrally coordinated magnetic semi- 
conductors has been put forward by the independent ef- 
fort of two . t e ams ., c onsisting of the present authors and 
co- workers.! 1 ^PIttI 13 ! In this theory, the ferromagnetic in- 
teraction between spins localized on the d shells of the 
magnetic ions is mediated by holes in the valence band. 
The free energy of the hole liquid is computed by di- 
agonalizing the 6x6 Luttinger Hamiltonian, which con- 
tains k ■ p, spin-orbit, and p-d exchange interactions, 
the latter taken into account in the molecular-field and 
virtual-crystal approximations. The influence of electro- 



static and spin disorder on magnetic properties is ne- 
glected, since they are not expected to have a qualita- 
tive impact on spin-polarized band electron thermody- 
namic properties. Hole-hole interactions can be taken 
into account in the spirit of Fermi liquid theory. The 
use of a mean-field approximation for the coupled band- 
electron and local-moment systems can be justified by 
the long-range character of the carrier-mediated spin- 
spin interaction, at least when the ratio of hole, to Mn 
density is sipall and the Fermi energy is large .lij It has 
been showr£3 that this model, with material parame- 
ters known from independent experiments, satisfactorily 
explains the magnitude of Tq, the temperature depen- 
dence of the spontaneous magnetization, the strength 
and strain dependence of the magnetic anisotropy, as 
well as the spectral dependence of the magnetic circular 
dichroism in (Ga,MpAAs. The hole-spin polarization has 
also been evaluated.E3 To address spin fluctuations a the- 
oretical description of ferromagnetism bevend the mean- 
field approximation has been de.ve]opedJl3 from which 
the magnon excitation spectrum, tBtB and the magnetic 
stiffness could be deduccd.113 

In this paper, we exploit this progres s , ia . the t heo- 
retical description_pf magnetic anisotropy£iJOli30 and 
magnetic stiffnessEj to address the domain structure in 
epitaxial layers of (Ga,Mn)As. The comparison between 
the computed and experimentally observed width of do- 
main stripes presented here constitutes an additional test 
of current theory. Moreover, the determined values for 
anisotropy energy, domain-wall energy and width may 
serve for—optimizing the design of, for instance, spin, 
injectionjij spin tunneling j 17 !! 18 ! or micromechanicalllj 
nanostructures of (Ga,Mn)As. 

For the calculation presented below we adopt band- 
structure parameters, elastic constants, and deformation 
potentials of GaA^-a set of values employed in our pre- 
vious work^O-BO for (Ga,Mn)As. The Mn spins are 
assumed to be in the d 5 configuration, so that S = 5/2 



and the Mn Lande. factor g = 2.0. For the p-d exchange 
energy we takeE3l!3 (3N = — 1.2 cV, which for the cation 
concentration of GaAs, N = 2.21 x 10 22 cm -3 , corre- 
sponds to J p d = —j3 = 0.054 eVnm 3 . The Fermi liquid 
parameter Ap = 1.2 enters the enhancement of Tq and 
of the p-d exchange splitting A = ApJ p dM/(g^B) rfthe 
valence band at magnetization M of the Mn spins. L3 

Another important parameter characterizing epitax- 
ial layers is the magnitude of biaxial strain. It de- 
pends on the layer thickness d and the difference be- 
tween the lattice parameters of the substrate and the 
layer, Aa = a s — a(x). In the case of (Ga,Mn)As that 
is obtained by low-temperatuxe MBE, films with d as 
large as 2 /zm are not relaxed.^ For such layers, the rel- 
evant components of the strain tensor assume the form 



magnetic theory is the width of the Bloch domain wall, 
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Cij are elastic constants. Since da(x)/dx = 0.032 nm 
the appropriately thin layer of (Ga,Mn)As deposited on 
GaAs or (Al,Ga)As is under compressive strain but the 
use of (Ga,In)As substrates can result in a tensile strain. 

It has been demonstrated by Ohno et ala that for 
films under compressive (tensile) strain the easy axis is in 
plane (perpendicular to film growth direction). Quanti- 
tatively, e xx = —0.2% for thg Gao.o965Mno.o3sAs film on 
GaAs, for which Ohno et al.a determined the magnitude 
of the magnetic field H a aligning magnetization along 
the hard axis. Tensile strain of e xx — 0.90% is expected 
for the 0.2 /im perpendicular film of Gao.957Mno.o43As 
on Gao.84Ino.i6As, a sample employed by Shono et alxA 
to examine the domain structure. Importantly, the the- 
ory referred to abovet^ILilli 2 ] reproduces correctly the 
strain dependence of magnetic anisotropy and, in par- 
ticular, explains thar4iiagnitude of H a .t3 Furthermore, 
it has been foundtlHl3 that for the relevant values of 
strain, hole concentrations p, and magnitudes of M, 
the energy density K u that characterizes uniaxial mag- 
netic anisotropy is greater than the corresponding cubic 
anisotropy terms K c and the energy density of the stray 
fields K d = fi M 2 /2. Thus, (Ga,Mn)As can be classified 
as a uniaxial ferromagnet. 

We discuss the d amain structure in terms of micro- 
magnetic theory.E3E3 Within this approach, a uniaxial 
ferromagnet is characterized by the anisotropy energy 
K Ul the magnetic stiffness A, and the saturation mag- 
netization M s . The procedure we adopt here consists of 
evaluating M(T), and thus A(T), in mean-field approx- 
imation. Then, K U (A) and ^4(A) are calculated. We 
believe that this procedure is well grounded at low tem- 
peratures. However fluctuation corrections will be more 
important at non-zero T, particularly in the critical re- 
gion near Tq. 

An important question arises whether the continuous- 
medium approximation underlying the micromagnetic 
theory is valid in diluted magnetic semiconductors, which 
contain a relatively low concentration of magnetic ions 
and even lower concentration of carriers. To address this 
question we note that the shortest length scale of micro- 



S w = ny/A/K u . 



(1) 



This length has to be compared to the mean distance 
between holes, r m = -y/3/(47rp). For concreteness, we 
consider Gao[»57Mno.o43As on Gao.84hio.i6 As for which 

10 20 cm" 3 we ob- 
91 K. The actual 



Tq = 80 K.B Assuming p = 3 
tain the mean-field value of Tq 
value of p is uncertain; the mean-field Tq would be 80 
K for p — 2.5 x 10 20 cm" 3 . Since the Fermi energy is 
greater than A at T — > (this corresponds to the "weak- 
coupling" or "RKKY" regime ,cJ) the redistribution of 
holes between the four valence subbands is only partial, 
and both A(A) and K U (A) are proportional to A 2 , ex- 
cept for large A, where the increase of K u with A is 
somewhat weaker. Hence, A/ K u and thus Sw are vir- 
tually independent of A, that is of T. For instance, for 
p = 3 x 10 20 cm" 3 , S w = 14.9 nm if T/Tq -» and 
decreases to 14.6 nm in the opposite limit T/Tq — > 1. 
Thus, the value determined for Sw ~ 15 nm is by more 
than a factor of ten longer than r m . We checked that 
Sw 3> r m in the whole relevant range of the hole concen- 
trations down to 10 20 cm" 3 . We conclude that the mi- 
cromagnetic theory in its standard continuous-medium 
form is suitable for modeling the domain structure in 
(Ga,Mn)As. 

In the case of perpendicular easy-axis films, a competi- 
tion between energies associated with the stray fields and 
the formation of the Bloch domain walls results in a sim- 
ple stripe domain-structure in the demagnetized thermal 
equilibrium state£3 Indeed, Shono et al. observed such a 
structure in Gao.957Mn .o43As/Gao.84l n o.i6As by means 
of a micro-Hall scanning probe.EU In particular, well de- 
fined stripes of width W increasing from 1.5 /iin at 9 
K to 2.5 /itm at 30 K were observed. Above 60 K, the 
stripes were less regular. Their width was evaluated to 
change from 3 to 6 ^m in the temperature range between 
65 and 77 K. The stripes were oriented along [110] at low 
temperatures but tended to lay along [100] abopjJ3p K. 

To interpret the experimental results we recalOcJ that 
for d ^> Sw and K u 3> Kd, the stripe width W is deter- 
mined by a solution of the transcendental equation, 



A c = (P 2 /7T 3 ) 



E 

n=l,3,5,.. 



n- 3 [l-(l + 27rn/P)e" 27r "/ p ]. (2) 



Here, P is the normalized stripe period P — 2W/d and 
the parameter A c describes the ratio of the Bloch domain- 
wall energy jw = ^\/AK u and the stray-field energy Kd, 



\ c = 4^/AKZ/(n M 2 d). 



(3) 



Figure 1 shows A c as a function of T/Tq computed for 
the film in question. Again the dependence on T is weak 
as both the denominator and numerator are to a good 
accuracy proportional to A 2 . 

We now evaluate the domain width W(T). The re- 
sults and a comparison to the experimental data of Shono 
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FIG. 1: Computed ratio of domain wall to stray-field en- 
ergy, A c (Eq. 5) as a function of reduced temperature for a 
Gai-zMn^As film with the easy axis along the growth direc- 
tion. 



et allEi are presented in Fig. 2. Furthermore, in order 
to establish the sensitivity of the theoretical results to 
the parameter values, we include results calculated for a 
value of A c 1.8 times larger as well. We see that the com- 
puted value for low temperatures, W = 1.1 fim, compares 
favorably with the experimental finding, W — 1.5 fim. 
The difference between the two may stem, in principle, 
from the uncertainty in the input parameters^- We note, 
however, that recent Monte Carlo simulations! 2 ^ demon- 
strates that disorder in positions of magnetic ions tends 
to enhance A, and thus W. At the same time, it is clear 
from Fig. @ that our model predicts much weaker tem- 
perature dependence of W than observed experimentally. 
It seems likely that the break in the experimental depen- 
dence W{T) around 60 K marks the beginning of the 
critical regime. In this regime, long-length-scale fluctu- 
ations in M, not accounted for in our theory, become 
important. These will have a larger effect on Kd, which 
is sensitive to fluctuations on scales shorter than W, than 
on the domain- wall energy 7^, which is sensitive only to 
fluctuations on length scales shorter than 5\y- 

Next, we comment on the crystallographic orientation 
of the domain stripes. The computed direction of the 
in-plane easy axis oscillates between [100] and [110] di- 
rections as aJiHLction of the hole concentration p and the 
splitting A.lilHi3 For the film in question we find that the 
lowest energy will have domain walls laying along [100] 
or equivalent directions, provided that p is greater than 
2.5 x 10 20 cm~ 3 . This stripe direction is in disagree- 
ment with the experimental resultscll described above. 
However, for slightly smaller values of p, the predicted 
behavior is in accord with the findings. In particular, 
for p = 2 x 10 20 cm -3 , the wall direction is expected to 
assume [110] or equivalent directions for large A — ► A s , 
that is at low temperatures, but to extend along [100] or 
equivalent directions at higher T, A < A s /2. 



FIG. 2: Temperature dependence of the width of do- 
main stripes as measured by Shono et alEi for the 
Gao.957Mno.o43As film with the easy axis along the growth 
direction (full squares). Computed domain width for the pa- 
rameter Ac depicted in Fig. hi is shown by the solid line. The 
dashed line is computed assuming that A c is by a factor of 1.8 
greater. 



It is well known that the domain structure of perpen- 
dicular easy-axis films exhibits an interesting evolution 
as a function of the magnetic field along the easy axis. 
When the field increases, stripes magnetized along the 
field grow, while those antiparallel to the field shrink. 
However, beyond a critical field value, cylindrical bubbles 
rather than stripes have lower energy. It would be inter- 
esting to search for such domains in perpendicular films 
of III-V magnetic semiconductors. At the same time, 
history-dependent metastable domain arrangements are 
expected to develop, for instance, a "froth" structure in 
the remanent state. The magnetic field H c \, at which do- 
mains vanish entirely under conditions of thermal equilib- 
rium increases with 1/A C , so that H c b — > M s for A c < 1. 
For the case under consideration, A c m 0.5 according to 
Fig. |], which corresponds to fJ, H c b « 0.1^ o M s « 5.5 
mT. This is consistent with the observed coercive force 
^ H C ~ 20 mT and square hysteresis for such a film.cd 
Actually, the fact that H c > H c i, implies the existence of 
a domain pinning mechanism. 

In view of the interest in magnetic nanocrystals, the 
length scale below which a ferromagnetic sample is in a 
single-domain state, is an interesting material parame- 
ter. Such a length scale is shape dependent. We consider 
square samples of ferromagnetic films with dimensions 
W x W x d. In thermodynamic equilibrium, the width 
Wsd below which the material is in a single-domain state 
is determined by d and l c = "fw /%Kd — A c d. For the per- 
pendicular film discussed above, d — 0.2 /im and l c = 0.1 
/im. For these values- according to numerical results of 
Hubert and Schafer,E3 Wsd — 1-2 /im. For such small 
single-domain particles, the celebrated Stoner-Wohlfarth 
theory predicts an abrupt switching of the magnetiza- 



tion direction in a magnetic field along the easy axis at 
Hf = H a = 2K u /(fi M s ). Thus, in this case the coercive 
force is equal to the anisotropy field aligning magnetiza- 
tion along the hard axis, /j, H c — fJ, H a = 670 mT for 
the material in question. It would be interesting to check 
experimentally the actual magnitude of H c in nanostruc- 
tures of III-V DMS. On the other hand, it is known al- 
ready that in the case of macroscopic films H c <C H a , re- 
flecting the inevitability of domain nucleation processes, 
associated-for instance-with space fluctuations of H a . 
The appearance of reverse domains followed by domain 
wall motion results in the complete reversal of magne- 
tization at H c in the range H a > H c > H c b. The wall 
motion begins when the field-induced torque on the wall 
magnetization overcomes the wall pinning force. 

In conclusion, our results imply that despite relatively 
small concentrations of magnetic ions and carriers, do- 
main properties of III-V DMS are, in many respects, 
similar to those of standard ferromagnets. In partic- 
ular, domain characteristic can be described in terms 
of micromagnetic theory. Such an approach, combined 
with our microscopic theory of hole-mediated ferromag- 
netism, predicts the direction and strength of magnetic 
anisotropy as well as characteristic dimensions of the do- 
mains correctly. We note that the values predicted for 



both anisotropy energies and domain-wall energies are 
qualitatively dependent on an accurate representation of 
the host semiconductor valence band. Further experi- 
mental studies on both macroscopic films and nanostruc- 
tures with fine spatial and time resolution will certainly 
improve our understanding of this novel ferromagnetic 
system, opening the doors wider for domain engineer- 
ing. Furthermore, the role of the magnetostriction and, 
in particular, its contribution to the strain tensor remains 
to be elucidated. On the theoretical side, it will be inter- 
esting to see to what extent intrinsic fluctuations in the 
Mn distribution, the distribution of other extrinsic de- 
fects, and the carrier density distributions they produce, 
account for domain nucleation and pinning fields. 
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